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ABSTRACT 

We present new sub-mm observations of Vega excess stars, and consistent modelling 
for all known Vega excess stars with sub-mm data. Our analysis uses dust grain models 
with realistic optical properties, with the aim of determining physical parameters of the 
unresolved disks from just their SEDs. For the resolved targets, we find that different 
objects require very different dust grain properties in order to simultaneously fit the 
image data and SED. Fomalhaut and Vega require solid dust grains, whilst HR4796 
and HD141569 can only be fitted using porous grains. The older stars tend to have 
less porous grains than younger stars, which may indicate that collisions in the disks 
have reprocessed the initially fiuffy grains into a more solid form, e Eri appears to 
be deficient in small dust grains compared to our best fitting model. This may show 
that it is important to include all the factors which cause the size distribution to 
depart from a simple power law for grains close to the radiation pressure blowout 
limit. Alternatively, this discrepancy may be due to some external infiuence on the 
disk (e.g. a planet). 

When the model is applied to the unresolved targets, an estimate of the disk 
size can be made. However, the large diversity in dust composition for the resolved 
disks means that we cannot make a reliable assumption as to the composition of the 
grains in an unresolved disk, and there is corresponding uncertainty in the disk size. 
In addition, the poor fit for e Eri shows that the model cannot always account for 
the SED even if the disk size is known. These two factors mean that it may not be 
possible to determine a disk's size without actually resolving it. 

Key words: circumstellar matter - dust. 



1 INTRODUCTION 

In 1983, the IRAS satellite detected a large infrared excess 
from Vega (a Lyr) during a routine calibration observation 
(Aumann et al. 1984). Soon, a similar excess was detected 
from three more stars (/3 Pic, Fomalhaut and e Eri). The 
excess could only be explained by a shell or ring of dust, 
heated to between 50 and 125 K by the central star. Quite 
by accident, IRAS had uncovered the first examples of solid 
material orbiting a normal main sequence star other that 
our own Sun. 

The four stars mentioned above have become the pro- 
totypes for the 'Vega excess' phenomenon (Backman & 
Paresce 1993; Lagrange, Backman & Artymowicz 2000). A 
systematic study of the IRAS all-sky survey shows that 
around 15 per cent of all main-sequence stars have an IR ex- 
cess due to circumstellar dust (Dominik & The Hjhvega Con- 
sortium 1997; Fajardo-Acosta et al. 1999). Coronographic 
and sub-mm imaging has shown the dust is typically in a 
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ring or disk geometry (Smith & Terrile 1984; Holland et al. 
1998) and that these systems have much in common with 
our own solar system, i.e. solid material in stable coplanar 
orbits distributed within a region roughly the size of the 
Kuiper belt. 

Though the Vega phenomenon is common, the poor res- 
olution of current far-IR and sub-mm instruments means 
that few Vega excess stars are close enough to be resolved. 
For the cases where the disks are resolved, detailed analysis 
using realistic models for the optical properties and size dis- 
tribution of the dust can be performed. This type of work has 
been attempted for /3 Pic (Li & Greenberg 1998a), HR4796 
(Augereau et al. 1999b; Li & Lunine 2003a) and Fomalhaut 
(Wyatt & Dent 2002). These models use the observed spa- 
tial distribution of dust, and fit the SED by varying the dust 
composition and size distribution. However, as these objects 
have been analysed on a case-by-case basis, different models 
and assumptions are used, which makes it somewhat diffi- 
cult to make direct comparisons. In addition, there are so 
few resolved targets that it is impossible to look for statisti- 
cal trends between disk size, age, stellar mass, etc. However, 
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there are many unresolved targets with good SED informa- 
tion, and new instruments such as SIRTF and SOFIA will 
produce many more such objects. A method which could re- 
trieve the physical parameters of a disk from just its SED 
would be invaluable. To accomplish this, the detailed prop- 
erties of the dust grains in Vega excess disks must be fully 
understood. 

We have tried to tackle these problems in two ways. 
Firstly, we have made sub-mm photometric observations of 
Vega excess stars, to help constrain their SEDs. Secondly, 
we have attempted to model all known Vega excess stars 
with sub-mm data, using a consistent model with realistic 
dust grains. This allows us to directly compare the results 
of modelling for resolved targets, and use the knowledge of 
dust grain size and composition gained to estimate the disk 
size for the unresolved targets. 



2 OBSERVATIONS 

Sub-mm observations were made using the SCUBA bolome- 
ter array on the James Clerk Maxwell Telescope. The ob- 
servations were made as part of several different observing 
programmes, and hence data was obtained in different ob- 
serving modes. In addition, we have reduced unpublished 
archive observations of suspected Vega-excess stars to en- 
large our sample. The observations are summarised in Ta- 
ble [H 

SCUBA can be used in two ways, either for photometry 
or for mapping. The distinction arises because the bolome- 
ters are spaced too widely to fully sample the image at the 
focal plane. Therefore in mapping mode, the moving sec- 
ondary mirror is used to add slight offsets to the telescope 
pointing, and hence fill in the gaps in the map. Whilst this 
is ideal for extended sources, for an observation of an un- 
resolved object a lot of time is spent where no bolometer 
can see the source. To prevent wasting observation time in 
this way photometry mode can be used, which does not at- 
tempt to make a fully sampled map. In this mode, the central 
bolometer is always on source, and hence the observation 
will be more sensitive. The disadvantage is that a fully sam- 
pled map is not produced, so if the object is resolved the 
measured flux will underestimate the true flux. This mode 
should therefore only be used if the target is known to be 
point-like. Although a full jiggle map is not used in this 
mode, a small 9-point map is used to prevent small pointing 
errors from affecting the measurement. 

For some of the observations an extended photometry 
mode was used, which used a 9-point map with a spacing of 
5 arcsec between each point, instead of the normal 2 arcsec 
spacing. The aim of this was to ensure that the true flux 
was measured even if the source was marginally resolved, 
whilst retaining most of the sensitivity advantage of pho- 
tometry mode. Also, some data was taken using a 12-point 
pattern, with four points at the centre and eight points ar- 
ranged around the central position, at a distance of 7 arcsec. 
Data from these extended photometry modes were reduced 
in exactly the same way as the conventional photometry. 



Object 


Date 


Wavelength 
(/im) 


Mode 


Hm 7206 


23 Aiip- 1 997 


450/850 


Map 




27 Aug 1997 


450/850 


Map 




10 Mar 2001 


450/850 


Ext. Phot. 


HD23362 


7 Jun 2000 


450/850 


Map 




2'^ Spd 1 997 


450/850 


Map 




29 Jan 2000 


450 /850 


Map 




21 Mar 2000 


450/850 


Map 




14 Aiiff 1997 


1350 


Phot 




13 Dec 1997 


450 /850 


Map 




14 Dec 1997 


450 /850 


Map 




1 1 Feb 1 998 


1350 


Phot 




29 Jan 2000 


450 /850 


Phot 


HD35187 


13 Dec 1997 


450/850 


Map 




14 Jan 1998 


450 /850 


Map 




11 Mar 2001 


450 /850 


Ext. Phot. 




12 Mar 2001 


450 /850 


Ext. Phot. 


HD48682 


9 Mar 2001 


450 /850 


Ext. Phot. 




10 Mar 2001 


450/850 


Ext. Phot. 




11 Mar 2001 


450/850 


Ext. Phot. 




12 Mar 2001 


450/850 


Ext. Phot. 


HD69830 


9 Mar 2001 


450/850 


Ext. Phot. 


HDSI ^1 ^ 

± ±J_^0 J- -LtJ 


21 Mar 2000 


450/850 


Map 


HD 109085 


21 Mar 2000 


450/850 


Map 




10 Mar 2001 


450/850 


Ext. Phot. 




11 Mar 2001 


450/850 


Ext. Phot. 




12 Mar 2001 


450/850 


Ext. Phot. 


HD121617 


19 Feb 1998 


1350 


Photometry 




11 Jun 2000 


450/850 


Photometry 


HD123160 


9 Dec 1997 


450/850 


Photometry 




19 Feb 1998 


450/850 


Photometry 




10 Jun 2000 


450/850 


Map 


HD128167 


21 Jan 2002 


450/850 


Ext. Phot. 


HD139664 


30 Nov 1997 


450 /850 


Map 




13 May 1998 


450/850 


Map 




25 Jun 1998 


450/850 


Map 




11 Jun 2000 


450/850 


Photometry 


HD141569 


18 Feb 1998 


1350 


Photometry 




19 Feb 1998 


1350 


Photometry 




4 May 1998 


1350 


Photometry 




13 Jan 1999 


1350 


Photometry 




10 Jun 2000 


450/850 


Photometry 




11 Jun 2000 


450/850 


Phot 


HD207129 


17 May 1998 


450/850 


Map & Phot 




28 May 1998 


450/850 


Map 



Table 1. New observations used in our analysis. Observations 
before 2000 were extracted from the JCMT archive. The extended 
photometry of HD128167 was observed with a 12 point pattern, 
all other extended photometry is 9 point. 

3 DATA REDUCTION 
3.1 Photometry data 

Photometry data were reduced using the SURF package. Af- 
ter combining the positive and negative beams, the data 
were flat-fielded and corrected for atmospheric extinction, 
estimated using a polynomial fit to the CSO 225 GHz opac- 
ity monitor and the relations in Archibald, Wagg & Jenness 
(2000). Residual sky emission was estimated using bolome- 
ters well separated from the source, and this signal was sub- 
tracted from the data. The sky bolometers were those in 
rings 2 and 3 (30 bolometers in total), and their median was 
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Object 


Other name 


Spectral type 


Distance 
pc 




Reference 


Estimated age 
Myr 


Reference 


HD17206 




F5/F6V 


13.97 










HD23362 




K2 


308.6 


7.9 X 10"'' 


1 






HD34282 




AO 


163.9 


0.39 


1 






HD34700 




GO 


> 180 


0.14 


2 






HD35187 




A2 


150 


0.14 


1 






HD38393 




F7V 


8.97 


5.3 X 10-6 


Best fit model 




8 


HD48682 




GOV 


16.4 


1.2 X 10"* 


Best fit model 










iVU V 


-LZ.UO 


2 3 X IQ— ^ 


Spline 


finn 9nnn 


1 n 






rtO V in. . . 
















F2V 


18 2 


4 X 10"^ 


opilllG 






HD121617 




AlV 


180" 


4.5 X lO"-^ 


2 






HD123160 




K5 


> 250 










HD128167 


cr Boo 


F2V 


15.5 


1.8 X 10-5 


Best fit model 


1700l"o° 


8 


HD139664 




F5IV-V 


17.5 


1.9 X 10-* 


Spline 


1120l«80 


8 


HD141569 




B9 


99.0 


8.3 X 10-^ 


2 


5±3 


7 


HD207129 




GOV 


15.6 


1.0 X 10-* 


Spline 




8 


HD22049 


e Eri 


K2V 


3.22 


1.1 X 10-* 


Best fit model 


730 ± 200 


10 


HD39060 


13 Pic 


A5V 


19.3 


3 X 10-3 


3 


20 ± 10 


6 


HD109573 


HR4796 


AOV 


67 


0.005 


4 


8±2 


9 


HD172167 


Vega 


AOV 


7.76 


2 X 10-5 


3 


354tf, 


5 


HD216956 


Fomalhaut 


A3V 


7.69 


8 X 10-5 


3 




5 



Reference list: 1 Sylvester et al. (1996); 2 Sylvester, Dunkin & Barlow (2001); 3 Backman & Paresce (1993); 4 Jura (1991); 
5 Song et al. (2001); 6 Barrado y Navascues et al. (1999); 7 Weinberger et al. (2000); 8 Lachaume et al. (1999); 9 Stauffer, 
Hartmann & Barrado y Navascues (1995); 10 Song et al. (2000). 

"No parallax was available, so the distance was estimated by assuming the absolute visual magnitude for a AlV star is 
+1.0 (Lang 1992, page 144), and using the measured visual magnitude of 7.298 from SIMBAD. 



Table 2. Parameters of stars used in our analysis. The top part of the table shows objects with new sub-mm data, the bottom part 
shows objects where we have analysed published data. Spectral type and distance are from SIMBAD. 



taken as the average sky signal. The signal from the central 
bolometer was then despiked using a 4 sigma cut. 

Calibration sources were observed and reduced in ex- 
actly the same way as the targets. If a planet was used as a 
calibrator, the fluxes package was used to estimate the flux 
and a correction was made for the loss of flux because the 
object was resolved. For the observations taken with the ex- 
tended photometry mode, all of the calibration sources were 
point-like (i.e. secondary calibrators, not planets). Thus the 
flnal photometry is derived under the assumption that the 
emission is point-like and centred on the star. If a target 
is extended (i.e. marginally resolved) or is offset from the 
position of the star, then observations in photometry mode 
will underestimate the true flux, by an amount shown in 
Table 13 As expected, the extended photometry modes are 
less affected than the conventional photometry mode. 



Offset 
arcsec 


FWHM 
arcsec 


Normal 


Relative response 

Ext. 9 pt. Ext. 12 pt. 








1.00 


1.00 


1.00 





5 


0.89 


0.93 


0.92 





10 


0.68 


0.76 


0.75 


5 





0.72 


0.81 


0.79 


5 


5 


0.67 


0.76 


0.75 


5 


10 


0.54 


0.65 


0.63 


10 





0.27 


0.41 


0.40 


10 


5 


0.28 


0.41 


0.39 


10 


10 


0.28 


0.39 


0.37 



Table 3. Comparison of the three photometry modes, showing 
effect of an extended or offset Gaussian source. The table shows 
the ratio of the quoted flux density to the true flux density. The 
PSF FWHM is assumed to be 14 arcsec. 



3.2 Maps 

Map data was reduced using custom data reduction software 
written in IDL. This system is similar to that developed to 
reduce data for the 8-mJy extragalactic survey (Scott et al. 
2002). This method is preferred over the SURF software be- 
cause the residual sky emission is removed more effectively 
(a plane is fltted to the residual emission, rather than a con- 
stant), and the 'zero footprint' rebinning method produces 
a map where each pixel is an independent measurement, 
which makes subsequent analysis much simpler. Source flux 
was measured by fitting an elliptical Gaussian to the source. 



To ensure that the software was reliable, SURF was used to 
reduce several sources, and the results from the two reduc- 
tion systems were compared to ensure that any differences 
were well within the errors. 

Calibration was done by calculating the fiux conversion 
factor from the peak flux of a Gaussian fit to calibration 
sources. Using the peak signal rather than the integrated 
signal means that the measurement is sensitive to changes 
in the beam shape (Dunne & Eales 2001), but this approach 
was necessary as most of the calibration sources were sec- 
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ondary calibrators, and only the peak signal is well deter- 
mined for these sources. 



4 RESULTS 

The results of the sub-mm observations are shown in Table|l| 
For objects with multiple observations, the flux given is the 
weighted mean of all the data. Of the 15 objects observed, 7 
were not detected, 4 were detected at 850 fim and 450 fim, 
and 4 were detected at 850 pm only. 

Most of the objects have not been observed in the 
sub-mm before, but for objects with existing observations, 
the new data were checked for consistency with the old 
results. HD34282 and H[D35187 have both been observed 
with UKT14 (Sylvester et al. 1996): HD34282 had a flux of 
1.3 ± 0.3 Jy at 400 fim and 0.409 ± 0.027 Jy at 800 ^m, 
HD35187 had a flux of 0.115 ± 0.022 Jy at 800 ^lm. The 
discrepancy between these fluxes and the new data is small 
(2 sigma), however similar discrepancies have been found 
in other comparisons of UKT14 and SCUBA data (e.g. 
Sylvester et al. 2001). This may be due to calibration er- 
rors in the original data. 

For the purposes of this paper, the extended photome- 
try observations have been reduced in exactly the same way 
as conventional photometry, so the result of the data reduc- 
tion is just the total flux. However, the data also contains 
some spatial information, and more complicated reduction 
will give some indication of the size of the disk as well as its 
flux. This will be tackled in subsequent paper (Wyatt et al., 
in preparation). 

One of the advantages of SCUBA observations is that 
they can confirm the IRAS excess is associated with the star, 
and not simply a background source (i.e. a galaxy or galactic 
cirrus). Contamination by background sources is a signifi- 
cant problem, and follow up observations are important to 
confirm that the IRAS associations are real. Certainly, not 
all IRAS associations are real, as demonstrated by observa- 
tions of 55 Cue (Jayawardhana et al. 2002) and HD155826 
(Lisse et al. 2002). SCUBA is a useful tool for confirming 
associations because it gives source positions which are sig- 
nificantly more accurate than those from IRAS. IRAS detec- 
tions have a typical positional uncertainty of about 16 arc- 
sec in the cross-scan direction and 3 arcsec in the in-scan 
direction (Beichman et al. 1988). In map mode, SCUBA po- 
sitional uncertainties are about 2 arcsec for bright sources, 
limited only by the pointing accuracy of the dish. In pho- 
tometry mode, the positional accuracy is determined by the 
size of the beam (FWHM of 14 arcsec) so these observations 
are less useful for this purpose. From our sample, one object 
was rejected on the basis of it's sub-mm position: HD123160 
was a strong SCUBA detection, but the source seemed to be 
offset from the stellar position by about 10 arcsec. In addi- 
tion, this star is now thought to be a distant class III giant 
(Kalas et al. 2002), and there is evidence of several galax- 
ies in the optical DSS survey image. This indicates that the 
IRAS source is probably a background galaxy. 

SCUBA is also sensitive to background galaxies, and so 
there is a possibility that a background galaxy could fall into 
the beam of the observations and cause a spurious detection. 
The chances of this happening depend on the limiting flux of 
the observations (as there are many more faint galaxies than 



bright ones) and the size of the beam. The extended photom- 
etry mode is therefore more susceptible to this problem, as 
it is sensitive to sources over ~ 2 times more area than con- 
ventional photometry. The tentative detection of HD38393 
has the highest probability of being due to a background 
galaxy, as the source is extremely faint and the observations 
were taken in extended photometry mode. This possibility 
is discussed further in Section [6.81 

The mass of the disks was estimated from the 850 ^im 
flux, using a dust temperature determined from the peak 
of the dust emission SED, and this is shown in Table 01 
We set the mass absorption coefficient k = 1.7 cm^ g~^ at 
850 fim for consistency with previous work (e.g. Zuckerman 
& Becklin 1993; Holland et al. 1998; Sylvester et al. 2001). 
Detailed discussion of the value of k can be found in Pollack 
et al. (1994). It is important to note that this mass esti- 
mate cannot measure the mass held in grains larger than 
about 1 mm. This is a key point, as in a collisional cascade 
(see Section 1^121 1 the total mass of the disk is dominated by 
the contribution from large bodies, which contribute little 
to the 850 ^im emission. Our mass estimate only refiects the 
amount of mass held in small dust grains, and even then is 
model dependent. The dust masses of the new targets are 
comparable with those measured in previous studies, rang- 
ing from 0.04 lunar masses (HD38393) to 33 lunar masses 
(HD141569) for the Vega excess stars, and hundreds of lunar 
masses or more for the pre-main-sequence stars (HD34700, 
HD35187 and HD34282). 



5 SED MODELLING 

Our key objective is to determine the size of a Vega excess 
disk from just its SED. Because models have many free pa- 
rameters, we have tried to constrain as much as possible us- 
ing our understanding of the physical processes taking place 
in the disks, in particular by setting the size distribution 
from theoretical arguments. 

The model used here is a modified version of the model 
developed to fit the SED and sub-mm image of Fomalhaut 
(Wyatt & Dent 2002), and a more detailed description can 
be found there. The model is based on a collisional cascade, 
where small dust grains are continuously created by colli- 
sions between larger bodies. The only significant difference 
between the model used here and that used in Wyatt & Dent 
(2002) is that we use a simpler assumption about the spatial 
distribution of the dust (see Section f5.4l for details). 

For this model to be accurate, the disk must be optically 
thin to radiation from the central star. If this is not the case, 
then the inner parts of the disk will shadow the outer parts 
from starlight, and a radiative transfer code would be needed 
to properly model the system. We have therefore excluded 
all the objects where Ldust/istar > 0.01. This ensures that 
disks are optically thin, as long as the opening angle is more 
than 2°. Measured opening angles are generally larger than 
this, e.g. 7° for 13 Pic (Heap et al. 2000). HD34282, HD34700, 
and HD35187 all have Ldust/^star > 0.1, so these objects 
have not been included in our analysis. 
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Object 


Flux (Jy) 


Tdust 


Dust mass 




450/im 


850Atm 


K 


.^^moon 


HD17206 


< 0.1789 


< 0.0078 




< 0.33" 


HD23362 


< 0.17 


< 0.012 




< 248° 


HD34282 


1.925 ± 0.29 


0.384 ±0.023 




3 X lO'* 


HD34700 


0.218 ± 0.037 


0.0407 ± 0.0024 




> 286° 


HD35187 


0.268 ± 0.047 


0.061 ± 0.005 




298° 


HD38393 


< 0.0428 


0.0024 ± 0.001 


82 


0.04 


HD48682 


< 0.0368 


0.0055 ±0.0011 


99 


0.3 


HD69830 


< 2.0248 


< 0.1005 




< 3.5° 


HD81515 


< 0.26 


< 0.02 




< 49.7° 


HD109085 


< 0.0381 


0.0075 ± 0.0012 


85 


0.6 


HD121617 


< 0.36 


< 0.012 




85 


HD123160'= 


0.152 ± 0.046 


0.021 ±0.007 






HD139664 


< 0.39 


< 0.015 




< 1.0° 


HD141569 


0.066 ± 0.023 


0.014 ±0.002 


90 


33 


HD207129 


< 1.7 


< 0.018 




< 0.95° 


a Boo 




0.0062 ±0.0017 


62 


0.5 


e Eri 


0.225 ± 0.010 


0.040 ± 0.0015 


85 


0.1 


Pic 




0.104 ±0.010 


103 


7.8 


HR4796 


0.180 ±0.15 


0.0191 ±0.0034 


99 


19 


Vega 




0.0457 ± 0.0054 


68 


0.8 


Fomalhaut 


0.595 ± 0.035 


0.097 ±0.005 


75 


1.7 



"Calculated assuming a dust temperature of 100 K. 
''Taken from Pietu, Dutrey & Kahane (2003). 

'^Sub-mm emission is offset by 9 arcsec from the star, and is thought 
to be due to a galaxy (see Section HI . 

Table 4. Photometry results and estimated dust masses for our sample. The top part of the table shows new data, the bottom part shows 
existing data. The dust temperature is estimated from the wavelength of maximum emission for the best fitting SED model. HD34700 
only has a lower limit on its distance, and hence we calculate a lower limit on its dust mass. These mass estimates are consistent with 
existing estimates made for Fomalhaut, Vega, and /3 Pic in Holland et al. (1998) and for e Eri in Greaves et al. (1998). 



Object 


12 /xm 


Flux (Jy) 
25 /jm 60 ^m 


100 Atm 


Fomalhaut 


12.1 


± 1.2 


3.9 ±0.3 


10.0 ±0.6 


10.7 ±0.6 


HR4796 






4.2 ±0.4 


7.1 ±0.9 


3.6 ±0.4 


Vega 


27.8 


± 1.4 


9.3 ±0.5 


9.8 ±0.4 


7.6 ±0.4 


t Eri 


6.6 ±0.3 


2.18 ± 0.11 


1.57 ±0.08 


1.99 ±0.22 


/3 Pic 






10.6 ±0.4 


19.1 ±0.8 


10.4 ±0.4 


HD38393 


3.00 


±0.15 


0.69 ± 0.04 


0.19 ±0.03 


< 0.5 


HD48682 


1.00 


±0.06 


0.29 ± 0.03 


0.36 ±0.05 


< 1.0 


HD141569 






2.29 ± 0.05 


5.3 ±0.8 


4.1 ± 0.5 


(J Boo 


1.10 


±0.05 


0.32 ± 0.02 


0.13 ±0.04 


< 0.5 


HD207129 


0.60 


±0.03 


0.15 ± 0.02 


0.31 ±0.04 


< 0.24 



Table 6. Colour corrected IRAS fluxes, based on the best fitting model. Upper limits are three sigma, and are taken directly from the 
IRAS faint source catalogue. HD109085 has been excluded because the best fitting model is unlikely to be representative of the true 
spectral shape at any of the IRAS wavelengths. 

5.1 Photometric data 5.2 Dust grain properties 

Our model of the dust grains is based on the assumption that 
the grains in protoplanetary disks are aggregates of interstel- 
Mid and far IR photometry data is taken from the IRAS lar dust grains, with an additional ice component frozen into 
faint source catalogue (Moshir et al. 1990). Where available, them as they grow. This type of model was originally devel- 
this is supplemented by ground based mid-IR data from the oped to simulate cometary grains (Greenberg 1982, 1998), 
literature. Also, published ISO photometry has been added but it has also been successfully applied to the dust grains 
to help constrain the far-IR part of the SEDs (in particular in accretion disks (Pollack et al. 1994) , and to the debris 
data from Habing et al. 2001). disks of (3 Pic (Li & Greenberg 1998a), HR4796 (Augereau 
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Object p gicc D^^i^ Reduced R 

fim AU 



Fomalhaut 


0.0 




9.6 


4.2 


150 




0.5 


0.0 


18.6 


14.1 


150 




0.9 


0.0 


91.0 


23.2 


150 




0.5 


1.0 


12.5 


9.3 


150 


HR4796 


0.0 




10.9 


16.8 


70 




0.5 


0.0 


21.3 


1.8 


70 




0.9 


0.0 


105.2 


3.2 


70 




0.5 


1.0 
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84 ±4 




0.9 


0.0 
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3.8 


43.8 ± 1.7 
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0.0 




4.2 


4.0 


180 ± 40 


HD38393 


0.0 




2.7 


0.4 


200 ± 50 


HD48682 


0.0 




2.3 


7.0 


71 ± 15"^ 




0.0 




2.3 


3.1 


110 ± 21'' 


CT Boo 


0.0 




3.3 


2.0 


320 ± 90 


HD207129 


0.0 




1.6 


3.7 


260 ± 50 



"Assuming secondary does not contribute to IRAS fluxes 
''Assuming secondary does contribute to IRAS fiuxes 



Table 5. Results of the SED fitting, where p is the porosity of the 
grains, q^^^ is the fraction of this porosity filled with ice, and -Dmin 
is the smallest dust grains in the system. -Dmin was determined 
by the blowout limit due to radiation pressure, except for e Eri 
(see Section 16.31 . Where no error is quoted on the disk radius, 
the size was constrained from existing images and was not a free 
parameter of the fits. 

et al. 1999b; Li & Lunine 2003a), HD141569 (Li & Lunine 
2003b) and HD207129 (Jourdain de Muizon et al. 1999). 

For the composition of the interstellar grains, we use 
the core- mantle model developed in Li & Greenberg (1997), 
which successfully fits the observed interstellar extinction 
and polarisation. This model has a silicate core surrounded 
by a UV-processed organic refractory mantle. We force the 
silicate/organic-refractory volume ratio to be 1:2, as inferred 
for cometary grains (Greenberg 1998). The grains are as- 
sumed to be porous, with the porosity being a free param- 
eter of the model. For dust grains cooler than ~ 120 K, the 
gaps due to porosity may be filled with either vacuum, wa- 
ter ice or a mixture of the two; for dust grains warmer than 
~ 120 K most ices would sublime, so the gaps must be filled 
with vacuum. The free parameters in the dust composition 



model are therefore the porosity p of the grains and the ice 
fraction gicc, with the constraint that for warm grains the 
ice fraction must be zero. 

Optical constants for the materials are taken from Li 
& Greenberg (1997, 1998a). The optical constants for the 
composite dust grain material are then calculated using 
Maxwell-Garnet effective medium theory. Then, the absorp- 
tion efficiencies for the dust grains are calculated using Mie 
theory, Rayleigh-Gans theory, or geometric optics, depend- 
ing on the size of the dust grain. Details of Maxwell-Garnet 
effective medium theory, and the methods used to calculate 
absorption efficiencies are discussed in Bohren & Huffman 
(1983). 

5.3 Size distribution 

The size distribution that results from an infinite coUisional 
cascade is (Dohnanyi 1969): 

n{D) oc D^-^" (1) 

where n{D)dD is the number of planetesimals of size be- 
tween D and D + dD, and q — 11/6. However, this distribu- 
tion will only hold for large dust grains, as radiation pres- 
sure will blow out small grains (e.g. < 1 nm) very rapidly, 
and Poynting-Robertson drag will cause intermediate sized 
grains (e.g. < 10 /im) to spiral inwards toward the star. 
Which of these two effects is dominant depends on the op- 
tical depth of the dust disk, but for all but the most ten- 
uous disks the coUisional timescale is shorter than the PR 
timescale, so PR drag can be safely ignored (Wyatt et al. 
1999). Given this, we assume that the above size distribu- 
tion holds down to the radiation blowout limit, when /3 > 0.5 
(where /3 = -FVad /-Fgrav ) , at which point there is a sharp cut- 
off. 



5.4 Spatial distribution 

For the spatial distribution of the dust grains, we assume a 
infinitely thin, flat ring, with the ring's radius R as n free 
parameter. This approach is motivated by the sub-mm im- 
ages of Fomalhaut and e Eri, which both show a narrow, well 
defined ring of emitting material. This spatial distribution is 
different from that used in Wyatt & Dent (2002), as here no 
attempt is made to account for the eccentricities of the dust 
grains, which broadens the ring by about 50 AU. The advan- 
tage of the simpler model is that it reduces the number of 
free parameters, but should still give a reasonable estimate 
of the overall extent of the disk. The Fomalhaut results from 
the two versions of the model are almost identical (see Sec- 
tion which gives us confidence that this simplification 
will not affect our conclusions. 



5.5 Stellar parameters 

Accurate stellar parameters are needed for several reasons. 
Firstly, the photospheric ffux must be estimated in each 
band, so that the excess ffux due to dust can be determined. 
This is particularly important in the IRAS 12 and 25 /im 
bands, where the photosphere can be the dominant source of 
ffux. Secondly, the temperature of the dust depends on the 



Suh-mm observations and modelling of Vega type stars 7 



luminosity of the star, and on the shape of the star's spec- 
trum in the region where most power is emitted (i.e. UV — 
optical — near-IR) . The luminosity and stellar spectrum also 
determine the magnitude of the radiation forces which act 
on the dust grains. Finally, an estimate of the mass of the 
star is needed to calculate the gravitational forces. 

To determine the stellar parameters, we first obtained 
the spectral type of each star from SIMBAD. This was used to 
estimate the star's effective temperature, using the calibra- 
tion provided in Gray & Corbally (1994). We then extracted 
a Kurucz model atmosphere with this effective temperature, 
assuming solar metallicity and surface gravity log((?) = 4.3 
appropriate for dwarf stars. The model atmosphere was then 
normalised to fit the K band magnitude. K band was pre- 
ferred over optical photometry because it is less susceptible 
to the effects of reddening. The luminosity of the star was 
then calculated by integrating the normalised model atmo- 
sphere. Finally, the mass of the star was estimated using 
the observed spectral type and the table provided by Lang 
(1992, page 132). 

5.6 Modelling method 

A model SED is specified by the grain composition (porosity 
p, ice fraction qicc), radius of the dust ring (R), and the total 
dust luminosity. In the cases where the size of the dust ring 
is known, the ring radius is fixed to the observed size; for the 
unresolved objects the ring radius is left as a free parameter. 
The approach we take is to choose a dust composition, and 
then find the best fit to the data by varying the remaining 
free parameters. 

The emission spectrum of an individual dust grain de- 
pends on its temperature and on the emission efficiency as 
a function of wavelength. The SED of a disk is simply the 
combined emission of all the dust grains. For a given dust 
composition model, the absorption efficiency is calculated as 
a function of grain size. Given this, the temperature of the 
dust grains as a function of grain size and distance from the 
star can be determined. Then, the total fiux as a function 
of wavelength is found by integrating over all grain sizes, 
weighting with the size distribution given in Equation 

The model spectrum is then converted into an SED that 
can be directly compared to the observed points to obtain 
a value. To do this, the spectral response of each filter- 
instrument combination is multiplied by the model spec- 
trum, to estimate the broad-band fiux that would be recov- 
ered from a real observation of the model spectrum. This 
step is particularly important for the IRAS data points, 
which have a very wide band-pass. By converting the model 
spectrum into broad-band fluxes, we avoid having to make 
colour-correction to the original data, which would other- 
wise require detailed knowledge of the true spectrum of the 
source. However, when showing the result of the fit, a colour- 
correction must be applied so that the data points can be 
plotted in real units. For this purpose, the best fitting model 
is used to estimate the colour correction for each point. If the 
best model is more than 3 sigma away from a data point, 
then no colour-correction is applied to that point, as the 
model is not likely to give a good estimate of the true colour- 
correction. 

For some objects it is impossible to fit all the data points 
simultaneously, because there is dust at a range of distances 
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Figure 1. SED of Fomalhaut and model fits with the ring radius 
set at 150 AU. Asterisks show optical and near-lR photometry, 
IRAS, ISO and SCUBA data are shown with error bars. Colour 
corrections for the IB. AS data were estimated from the p = 
model. 

from the star. In these cases, only the long wavelength points 
are fitted, in order to obtain an estimate of the overall extent 
of the ring. 



6 MODELLING RESULTS 
6.1 Fomalhaut 

The size of Fomalhaut 's dust disk is known from the sub-mm 
images made with SCUBA (Holland et al. 2003). Because 
the ring is inclined and the images have low resolution the 
size of the ring cannot be measured directly from the images, 
but the best estimate of the true size is 150 AU (Wyatt & 
Dent 2002). The SED data for Fomalhaut comes from IRAS, 
SCUBA and also ISO. The ISO fiuxes are from Walker et 
al. (in preparation), and are listed in Wyatt & Dent (2002). 

Our modelling (Figure^ shows a good fit for solid dust 
grains (i.e. p = 0.0), but rules out porous grains as this 
would make the dust hotter than is observed. The p = 0.5, 
Qicc = 1.0 model is also ruled out, as it fails to fit the mid-IR 
part of the SED. 

The reduced for the best fitting model is 4.2, which 
indicates that the model is not consistent with the data. 
However, by far the dominant contribution to comes from 
the ISO data, which show considerable scatter from any 
smooth model. The ISO filters all have broad band passes 
(A A/A ~ 0.5), so the apparent scatter in the ISO points 
is unlikely to be due to structure in the spectrum. A more 
likely explanation is that uncertainties on the ISO points 
have been underestimated, and this would explain the high 
reduced x^- 

Fomalhaut has already been modelled in detail by Wy- 
att & Dent (2002), and the results here are almost identical 
to those results. As the modelling here uses the same data 
and dust grain model, it is hardly surprising that we re- 
cover the same results. However, this does confirm that the 
simpler spatial distribution used here does not significantly 
affect the results of the modelling. 
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Figure 2. SED of HR4796 and model fits with the ring radius set 
at 70 AU. Data where the wavelength is shorter than 18 fim has 
been excluded from the fit, as the excess here is probably due to 
a hot dust component close to the star (see Section ^6.21 . Upper 
limits are 3 sigma. Colour corrections for the IRAS data were 
estimated from the p = 0.5, 5ice = model. 



6.2 HR4796 

HR4796 is resolved in both mid-IR images of dust emission 
(Jayawardhana et al. 1998; Koerner et al. 1998) and scat- 
tered light images in the near-IR (Schneider et al. 1999). 
The best estimate of the size of the ring comes from mid-IR 
images, which suggest that the radius is 70 AU. The star is 
younger than most of the other stars in the sample, with a 
age of 8 ± 2 Myr (Stauffer et al. 1995). 

The modelling results show that solid dust grains are 
unable to fit the observed SED and disk size (Figure I^J. 
However, using a porosity of p = 0.5 provides a good fit to 
the SED. The icy model {p — 0.5, gicc ~ 1.0) is a worse 
fit than the solid grain model, so can be rejected. Around 
10 fim, there is an additional excess which cannot be ac- 
counted for by any model. This is likely to be produced by 
a hot dust component close to the star. Data in this part of 
the spectrum have therefore been excluded from the fit, and 
do not contribute to the total quoted in Table |S] 

Detailed modelling of this object has been done previ- 
ously by Augereau et al. (1999b), and our results here are 
in good agreement with this work. Augereau et al. made a 
much more thorough exploration of parameter space, but 
their overall conclusion was that two dust populations are 
required: a cool component at 70 AU, and a hotter compo- 
nent at about 10 AU. Their best fitting model for the cool 
dust has amorphous grains with a porosity of p = 0.6 and 
an ice fraction ice gico = 0.03. In their analysis, the ring's 
diameter was constrained to be 70 AU, but the grain poros- 
ity, ice fraction and the width of the dust ring were all left 
as free parameters. As for Fomalhaut, our results confirm 
that the simpler model used here is able to produce similar 
results to a more detailed model. 

Modelling of HR4796 has also been done by Li & Lunine 
(2003a), and they were able to fit the entire SED without 
having a hot dust component, in contrast to the results of 
Augereau et al. (1999b) and the results presented here. This 



Figure 3. Graph showing how /3 (= i^rad/^grav) depends on 
dust grain size, for dust orbiting e Eri. /3 = 0.5 is shown by the 
horizontal line; any dust grain with /3 > 0.5 will leave the system 
on a hyperbolic trajectory. 

discrepancy is probably because Li and Lunine's model con- 
tains a large number of grains which are much smaller than 
the blowout diameter, and these small grains produce a hot 
component in the SED even though they are distant from 
the central star. 



6.3 e Eri 

e Eri is a particularly interesting object, as it is the only 
cool star (K2V) with a resolved dust disk, all other resolved 
disks being around A stars. SCUBA images of e Eri show 
a distinct ring, which is approximately face on and has a 
radius of 60 AU (Greaves et al. 1998). 

In modelling e Eri's dust disk, a problem arises due to 
the low luminosity of the star. Our model for the dust size 
distribution assumes that radiation pressure will remove the 
smallest dust grains. However, in the case of e Eri the radia- 
tion pressure is not sufficient to blow dust grains out of the 
system when porous dust grains are used (Figure EJ. In this 
situation, it is not clear what limits the size of the smallest 
dust grains. As a first step, we made the ad-hoc assumption 
that when porous grains were used, the small size cutoff is 
set at the grain size where /3 is greatest, i.e. where the radi- 
ation force is most significant compared to the gravitational 
force. Modelling using this assumption produced a very in- 
teresting result: all of the models predict too much flux in 
the mid-IR, and too little in the sub-mm (Figure 0J . This 
means that our models have too many small (and therefore 
hot) dust grains. This is an unexpected result, as there is 
no obvious mechanism for removing the small grains. By 
varying the minimum size cutoff explicitly, an excellent fit is 
possible if the cutoff size is set to 3.5 /im and the p = model 
is used, as shown in Figure |3 It is also possible to obtain a 
good fit using porous grains, but only if the minimum size 
cutoff is set to around 300 fim (Figure (SJ. This model does 
not predict the observed excess at 25 fim, but this flux could 
come from an additional hot component close to the star. 
Modifying the slope of the size distribution can improve the 
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fit (Figure |7||, but cannot make the model consistent with 
the data. 

Sohd dust grains with a size of 3.5 fim have a /3 of about 
0.15, much less than is needed to remove the grains by radi- 
ation pressure. It therefore seems that the disk around e Eri 
contains few grains smaller than about 3 /xm, and radiation 
pressure from the star cannot explain this deficiency. Possi- 
ble causes for this discrepancy are discussed in Section |7| 

Recently, detailed modelling of e Eri has been pre- 
sented by Li, Lunine & Bendo (2003) who modelled the SED 
and SCUBA image using a porous dust model {p = 0.9, 
5icc = 0.0), and were able to get an excellent fit. With 
this in mind, our results seem somewhat surprising, as they 
favour solid rather than porous grains, and fail to produce 
an acceptable fit at all unless grains with diameter smaller 
than 3.5 ^m are excluded. The discrepancy arises because 
of several differences in the assumptions made within the 
two models. The most important difference is that Li et al. 
(2003) set a maximum grain size cutoff at a diameter of 2 cm, 
whereas in our modelling the maximum grain size is set to 
be 10 m. If the collisional cascade model is true, then the 
real size distribution will extend up to very large bodies (i.e. 
kilometre scale and bigger), but the limit of 10 m is chosen 
because bodies of this size and larger make a negligible con- 
tribution to the sub-mm flux (i.e. the longest wavelength 
data included in our fits), and so can be safely excluded 
from the model. This is not true of 2 cm bodies, as grains 
of around this size can make a very significant contribution 
the the flux at 850 fj,m. Similarly, Li et al. (2003) set an a 
priori limit on the minimum grain diameter at 2 fim; our 
flnding that we must remove grains smaller than 3.5 fim is 
therefore in reasonable agreement with this (the remaining 
discrepancy is due to the fact that Li et al. (2003) allow the 
slope of the size distribution to be a free parameter, whereas 
in the models where we modified 73min, q was fixed to be 
11/6). Finally, Li et al. (2003) have a different spatial dust 
distribution, with a peak at 55 AU and FWHM of 30 AU, 
compared with the infinitely thin ring at 60 AU used in our 
model. The broad spatial distribution has not been consid- 
ered here as it would make our SED fit worse. Because Li 
et al. (2003) set tight limits on the size distribution, the dust 
SED becomes quite sharply peaked, and a broad spatial dis- 
tribution is needed to make the dust SED wide enough to fit 
all the points. Conversely, in our model the broad size dis- 
tribution means that the model SED is already too wide, so 
a broad spatial distribution would only make things worse. 
Given the low signal to noise of the SCUBA image, a direct 
measurement of the width of the dust ring is not possible. 

Dent et al. (2000) also modeled the e Eri data, and were 
able to fit the SED using a modified blackbody. The fitted 
parameters indicated that the disk must contain large grains 
(~ 30 fJ,m), but this type of modelling does not constrain the 
composition or porosity of the dust grains. 

6.4 Vega 

The original estimate for the size of Vega was made by Au- 
mann et al. (1984), who analysed special pointed observa- 
tions made with IRAS after the discovery of a far-IR ex- 
cess. The estimated diameter was 23 arcsec at 60 /xm (cor- 
responding to a radius of 90 AU). Subsequently, this data 
was reanalysed, giving a larger diameter of 35 ± 5 arcsec 
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Figure 4. SED of e Eri, with the ring radius set at 60 AU. Colour 
corrections for the IRAS data as plotted here were estimated from 
the p = model. 
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Figure 5. SED of e Eri, with the ring radius set at 60 AU. Models 
show the effect of changing the minimum size cutoff when using 
solid grains (p = 0). Colour corrections for the IRAS data were 



estimated from the D„ 



■ 3.5 /xm model. 



(135 AU) (van der Bliek, Prusti & Waters 1994, and refer- 
ences therein). This discrepancy is attributed to the fact the 
Aumann et al. did not account for the photospheric emis- 
sion at 60 ^m, and hence underestimated the extent of the 
excess emission. However, sub-mm observations made with 
SCUBA on the JCMT indicate a diameter of 24 ± 3 arc- 
sec (Holland et al. 1998), similar to the original estimate 
from the IRAS data. At millimetre wavelengths, aperture 
synthesis imaging has been done using the Plateau de Bure 
interferometer (Wilner et al. 2002) and Owens Valley Ra- 
dio Observatory (Koerner, Sargent & Ostroff 2001). Both 
of these observations suggested that the Vega disk is very 
clumpy, and found bright clumps at a distance of 9 arcsec 
(Plateau de Bure) and 12 arcsec (Owens Valley) from the 
central star. This could indicate a disk radius as small as 70 
AU. 
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Figure 6. SED of e Eri, with the ring radius set at 60 AU. Models 
show the effect of changing the minimum size cutoff when using 
porous grains (p = 0.9). Colour corrections for the IRAS data 
were estimated from the -Dmin = 300 fim model. 
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Figure 7. SED of e Eri, with the ring radius set at 60 AU. Models 
show the effect of changing the slope of the size distribution (q). 
Colour corrections for the IRAS data were estimated from the 
q = 1.8 /im model. 



Given these observations, it is not obvious what the true 
size of the disk is, and so we were unable to fix this parame- 
ter in our model fits. Instead, we fitted different grain mod- 
els with disk size as a free parameter. With this approach 
we find that all of the dust grain models fit equally well 
(Figure ISJ, but predict very different disk sizes. The porous 
grain models both have a fitted radius much larger than any 
of the observations, with the p — 0.5, gicc = 0.0 model sug- 
gesting a radius of 191 ± 16 AU, and the p = 0.9, gico = 0.0 
model giving 237 ± 21 AU. These models can therefore be 
rejected. However, the solid grain model {p — 0.0) gives a 
radius of 120 ± 13 AU, and the icy p = 0.5, gicc = 1.0 model 
gives 92 ± 8 AU, both of which are compatible with the ob- 
servations. Without more detailed analysis of the structure 



Figure 8. SED of Vega and model fits to the data with ring 
radius as a free parameter. Colour corrections for the IRAS data 
as plotted here were estimated from the p = model. 
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Figure 9. SED of /3 Pic and model fits to the data with ring 
radius as a free parameter. Colour corrections for the IRAS data 
as plotted here were estimated from the p = model. 



observed in the disk, it is impossible to distinguish between 
these two options. 

There is a small additional excess at 25 fim ('^ 2 Jy), 
which none of the models account for. This could either 
be due to an additional warm dust component similar to 
that found for HR4796, or because the spatial distribution 
is broad, with dust at a range of distances from the star. 



6.5 13 Pic 

/3 Pic is the most extensively studied Vega excess star, with 
observations in scattered light, mid-IR, far-IR, sub-mm. The 
system is not a ring, but a disk with dust at a variety of 
distances from the star. The system has been modelled in 
detail by Li & Greenberg (1998a), and our simpler model is 
not adequate to deal with the extended spatial structure of 
the disk. However, we have included f3 Pic here in order to 
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show the consequences when our model is applied to objects 
where the dust is not contained within a ring. The ring mor- 
phology is one of the key assumptions of our model, so it is 
extremely important that we can determine whether or not 
this is true for unresolved objects. 

To fit the SED, we used the same approach taken for 
Vega, i.e. we left the ring radius as a free parameter, and fit- 
ted each dust composition to the SED. The modelling shows 
that it is not possible to fit the observed SED with any single 
model regardless of what dust composition is chosen (Fig- 
ure (nj. This is a helpful result, in that it shows that we are 
able to distinguish between a disk and a ring morphology 
purely from the SED. However, this modelling reveals lit- 
tle else about the composition and size distribution of dust 
grains in (3 Pic. 
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6.6 HD141569 

HD141569 has a resolved circumstellar disk, first directly 
detected in near-IR scattered light images (Augereau et al. 
1999a; Weinberger et al. 1999) , and subsequently resolved in 
mid-IR thermal emission (Fisher et al. 2000; Mouillet et al. 
2001; Marsh et al. 2002). Most recently, the disk has been 
observed with the ACS coronagraph (Clampin et al. 2002), 
revealing that the structure observed in the disk is probably 
caused by a tidal interaction with a bound binary system at 
a projected distance of around 1000 AU. The age of the sys- 
tem has been estimated to be 5 ± 3 Myr (Weinberger et al. 
2000), making it one of the youngest stars in our sample. 
However, despite its youth the disk appears to be in a col- 
lisional cascade, as modelling shows that it contains small 
dust grains which must leave the system on a very short 
timescale (Fisher et al. 2000). 

The mid-IR imaging shows that the overall size of the 
disk is around 100 AU, but the emission is produced at 
a range of distances from the star from 20 AU outwards 
(Marsh et al. 2002). This is also apparent from our SED 
modelling, as no model can simultaneously fit all the data 
(Figure EU- When only the long wavelength points are fit- 
ted (i.e. only the cooler, outer parts of the disk), the fitted 
radius is about 50 AU if solid (p = 0) or icy (p = 0.5, 
Qicc ~ 1.0) grains are used, but 114 AU for porous grains 
(p = 0.9, Qico = 0.0). This indicates that the thermal emis- 
sion is probably produced by porous grains, in agreement 
with the modelling done by Li & Lunine (2003b). 

6.7 HD109085 

HD109085 is clearly detected in the sub-mm, indicating that 
the IRAS association is real. However, no model can simul- 
taneously fit all of the data (Figure [TT|l . This indicates that 
there must be dust at a variety of distances from the star, as 
is the case for /3 Pic and HD 141569. For our fitted model, the 
IRAS 12 and 25 pm points have not been included, so the fit- 
ted radius refiects the size of the coolest parts of the disk, i.e. 
the disk's overall extent. This gives a result of 180 AU, which 
at a distance of 18.2 pc gives an angular radius of 9.9 arcsec. 
This is more extended than the SCUBA beam, so our flux 
estimate is probably smaller that the true flux by a factor 
of ~ 2 (Table 1^. Given the size and brightness of the disk, 
it should be possible to fully map HD109085 with SCUBA, 



Figure 10. SED of HD14159 with model fits for different grain 
compositions. The ring radius is a free parameter in each case. 
Colour corrections for the IRAS data as plotted here were esti- 
mated from the p = model. 
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Figure 11. SED of HD109085 with model fit to the data. Ring 
radius is a free parameter. Upper limits are 3 sigma. 

and such observations are currently being performed (Wyatt 
et al. in preparation). 

6.8 HD38393 

HD38393 is extremely faint in the sub-mm, and is only de- 
tected at the 2.4 sigma level. The modelling produces a 
good fit to the SED (Figure I12II . but with a disk radius 
of 200 ± 50 AU for solid dust grains. At the distance of 
HD38393 (9.0 pc), this corresponds to a radius of 22 arcsec. 
This would place most of the flux outside the SCUBA beam 
(even though the extended photometry mode was used), im- 
plying that the true flux at 850 nm is much higher than the 
measured flux. In fact, if this radius is correct we would only 
detect emission if the source is edge on. Revising our esti- 
mate of the 850 fim flux results in a still larger fitted radius 
(because a higher sub-mm flux indicates a cooler disk), and 
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Figure 12. SED of HD38393 with model fit to the data. The ISO 

60 /im point has been offset to the right to distinguish it from the 
IRAS 60 /im point. Upper limits are 3 sigma. 



Figure 13. SED of a Boo. The 750 60 /im point has been oflFset 
to the right to distinguish it from the IRAS 60 /im point. Upper 
limits are 3 sigma. 



so docs not resolve this problem. Nor does choosing a diflter- 
ent dust composition, as porous grains tend to be warmer 
than solid grains, and so also increase the estimated disk 
size. 

Given the low significance of the detection, it is possible 
that the 850 /xm detection is either spurious, or caused by a 
background galaxy. Using source counts of galaxies detected 
in blank fields (e.g. Figure 11 in Scott et al. 2002), we find 
there is a 6 per cent probability that a 2.4 mjy background 
source will fall into the beam of SCUBA iu extended pho- 
tometry mode (which is sensitive over a 300 arcsec^ area). 
Habing et al. (2001) also concluded that their 170 jim de- 
tection is likely to be due to a background galaxy, given the 
low flux and the galaxy rmmber counts at that wavelength. 
However, an alternative explanation is that the estimate of 
the disk size is wrong: the difficulty fitting the SED of e Eri 
may indicate that our model is generally unreliable for cool 
stars, and hence the true size of the disk could be smaller 
than 200 AU. This would mean that the disk is deficient in 
small grains compared with our model. 



6.9 a Boo 

Excess sub-mm emission from a Boo was detected at the 3.6 

sigma level, with a flux of 6.7 mJy. However, the best fitting 
model gives a disk radius of 320 ±90 AU, which corresponds 
to an angular radius of 20 arcsec, implying that most of 
the 850 /im flux would be outside the SCUBA beam even 
in extended photometry mode. In addition, the best fitting 
model predicts a much lower 850 /im flux than the observed 
value. As is the case for HD38393, there is no model which 
is consistent with both the SED and the spatial constraints 
given by the SCUBA beam size. The problem is unlikely 
to be caused by errors in the far-IR data: the IRAS and 
ISO results arc entirely consistent, and the ISO 170 /im 
observation has a beam radius of about 45 arcsec, so all the 
flux would be detected with ISO, even if the disk really is 
very large. 



Habing et al. (2001) considered the possibility that 

their ISO detection at 60 and 170 /im was duo to a back- 
ground galaxy. They concluded that the excess was probably 
real, but could not rule out contamination by a background 
source. A background galaxy with a flux of 6.7 mJy is un- 
likely to fall into the SCUBA beam by chance alone (1 per 
cent probability), but as a Boo was selected on the basis of 
a far-IR excess, contamination by a background source be- 
comes more likely. The alternative explanation is that a Boo 
is deficient in small grains, which would reconcile the mod- 
elling results with the SCUBA beam size. 



6.10 HD48682 

HD48682 is a quoted as visual binary in the Washing- 
ton double star catalogue with a separation of 34 arc- 
sec. However the two stars cannot be physically associated 
due to their differing proper motions, and it appears that 
HD48682B is in fact a background object. Given that the 
separation of the two components is similar to the IRAS 
beam size, it is not obvious whether the photosphere of the 
secondary will contribute to the IRAS fluxes or not. The 
IRAS beam is different at each wavelength, as the diffrac- 
tion limit is much smaller at 12 /im than it is at 100 /im. This 
means that the secondary photosphere may only contribute 
to the measured fluxes at longer wavelengths. In addition, it 
is very difficult to estimate the brightness of the secondary 
at IRAS wavelengths, because the star is most likely an M 
star and only optical photometry is available. A small uncer- 
tainty in the effective temperature therefore causes a large 
error in the estimated mid-IR flux. A K band image would 
resolve this problem, but unfortunately the secondary star 
is saturated on 2MASS images, so a dedicated observation 
would be required. Aumann & Probst (1991) made ground 
based observations of HD48682 at 10 /im, and concluded 
that the 12 /im excess comes from > 6 arcsec from the pri- 
mary and is consistent with the flux from the secondary, but 
that the secondary photosphere cannot account for the large 
60 /im excess. Our detection with SCUBA confirms that the 
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Figure 14. SED of HD48682, with model fit to the data. As- 
terisks show optical and near-lR photometry, IRAS and SCUBA 
data are shown with error bars. The photospheric flux from the 
primary, and the best fit model assuming that only the primary 
contributes to the IRAS fluxes are shown with dashed lines. The 
secondary is shown with a dotted line. Solid lines show the com- 
bined photospheric flux, and combined photospheric flux plus best 
fit model assuming that the secondary does contribute to the 
IRAS fluxes. Upper limits are 3cr. 

60 fj,m IRAS excess is associated with the primary and not 
the secondary, but the source of the 12 and 25 fim excess is 
not clear. 

Figure ITU shows the SED of HD48682, and the best fit 
model. The modelling has been done in two ways, with and 
without the secondary included in the photosphere subtrac- 
tion (assuming a spectral type MO). At 12 ^m, the IRAS flux 
is more than the photosphere of the primary, which could ei- 
ther be due to circumprimary emission at > 6 arcsec, or flux 
from the secondary. At 25 fj,Tn, the IRAS flux is more than 
the combined flux from primary and secondary, indicating 
that there is definitely emission by dust at this wavelength. 
There is a large excess at 60 fim, and the secondary pho- 
tosphere could only make a small contribution even if it is 
within the beam. The modelling results are better if the sec- 
ondary is included in the photosphere subtraction (reduced 
= 3.1, verses 7.0 if the secondary is not included). If the 
secondary is not included in the subtraction, the estimated 
radius is 71 ± 15 AU, whereas if the secondary is included, 
the estimated radius is 110 ± 21 AU. If the larger value is 
true, the projected radius would be 6.7 arcsec, so the disk 
could be marginally resolved with SCUBA mapping. The 
shape of the fitted model indicates that there would be only 
a small excess at 18 fj,m, so resolved ground based imaging 
would be difficult. 

6.11 HD207129 

HD207129 is not detected in our reduction of archived 
SCUBA measurements, but because there are ISO measure- 
ments at 170 a model fit is possible, and the SCUBA 
upper limit also helps to constrain this fit. The SED and 
model fit are shown in Figure \TEl The fitted size of the disk 
is 260 ± 50 AU, which corresponds to an angular radius of 



Figure 15. SED of HD207129. The ISO 60 fim point has been 
offset to the right to distinguish it from the IRAS 60 fira point. 
Upper limits are 3 sigma. The lower data point at 850 /ivn is the 
result if the source is unresolved, the upper point is the result if 
the source is resolved but with a radius smaller than 20 arcsec. 

17 arcsec. Since this is larger than the SCUBA beam a sig- 
nificant fraction of the flux would not have bean detected, 
and the quoted upper limit on the 850 flux may actu- 
ally be lower than the true flux. HD207129 was observed 
in both photometry mode and map mode, and the quoted 
upper limit of < 0.018 Jy is based on both sets of measure- 
ments. If just the mapping data is used, then an integrated 
flux over a 20 arcsec radius aperture can be calculated. This 
method gives a flux < 0.05 Jy (3 sigma). 



7 DISCUSSION 

The most striking result of our modelling is that there is a 
large diversity in dust grain compositions between different 
disks. A simple explanation for this result is that dust grain 
composition may vary with age. The modelling of HD141569 
and HR4796 (age 5 and 8 Myr respectively) indicates high 
porosity grains (p = 0.9), as does the previous modelling 
of P Pic (age 20 Myr) done by Li & Greenberg (1998a). 
However, solid grains are indicated for Fomalhaut and Vega 
(ages 150 and 350 Myr). It may be that the older stars have 
less porous grains because the collisional cascade has repro- 
cessed the initially fluffy grains into a more solid form. An 
alternative explanation is that the change in porosity is due 
to the effects of stellar radiation on the dust, as suggested 
to explain the differences between the porosity of long and 
short period comets (Mukai & Fechtig 1983; Li & Greenberg 
1998b). Clearly, a larger sample of resolved disks is required 
to confirm this trend. 

The failure of our basic model to fit the SED of e Eri is 
extremely interesting, and it's important to determine what 
is causing the poor flt. There are essentially two possibilities; 
either there is an unknown physical effect which is modifying 
the size distribution in the particular case of e Eri, or else 
the size distribution model we have assumed is intrinsically 
wrong. 

e Eri is known to have a planet (Hatzes et al. 2000), 
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with a semi-major axis of 3.4 AU and a mass of around 
2 Jupiter masses (assuming an inclination of 25° as measured 
by Greaves et al., 1998). This planet in unlikely to affect the 
dust size distribution at 60 AU, because of its small orbital 
radius. However, a planet with a large semi-major axis has 
been inferred from the existence of clumps within the dust 
disk (Quillen & Thorndike 2002). If this planet is real, then 
dynamical effects axe clearly very important in shaping the 
disk, and this would undoubtedly affect the size distribution 
of the dust (Wyatt 2003, submitted). 

The alternative explanation is that the discrepancy is 
not due to an external influence, but that a more accurate 
treatment of the dust grain size distribution in a coUisional 
cascade is needed to account for the SED. More detailed size 
distribution models allow the grain strength to vary with 
grain size, and account for the the effects of removing small 
grains on the overall size distribution (Campo Bagatin et al. 
1994; Durda & Dermott 1997). For Pic, extensive rmmer- 
ical models of this type have been produced, both for the 
inner disk (Thcbault, Augereau & Beust, submitted) and for 
the outer disk (Krivov, Mann & Krivova 2000) . Both of these 
studies show significant deviation from the theoretical size 
distribution used in this paper. This suggests that a more 
rigorous treatment of the coUisional cascade in e Eri would 
also give a significantly different size distribution. However, 
it is not obvious that this would produce a better fit to 
the SED (i.e. that the new size distribution would have 
few grains smaller than 3 fiin). In fact, according to Campo 
Bagatin et al. (1994) we should expect an enhancement of 
grains just above the blowout radius as they are less likely 
to be destroyed by collisions, which would make the prob- 
lem worse. Conversely, the modelling of /3 Pic's outer disk 
by Krivov et al. (2000) suggests a reduction in the number 
of grains just above the blowout limit, because they are de- 
stroyed by small, fast moving grains as they are ejected from 
the system by radiation pressure. However, as the stellar lu- 
minosity is much lower for e Eri, this mechanism is likely to 
be less important. Only detailed modelling specific to e Eri's 
dust disk can answer this question. 

The new sub-mm data significantly expands the sample 
of true Vega excess stars detected at this wavelength. Re- 
gardless of the dust properties, we are able to distinguish 
between a disk morphology and a ring morphology. Out of 
the five new objects, only HD109085 is disk like, indicating 
that this type is less common than the ring like structure 
observed around Fomalhaut and e Eri. 

One of the main aims of this work was to determine the 
size of an unresolved disk from just its SED. However, the 
modelling of resolved targets has shown that this may not 
be possible. The large diversity in dust properties for the 
resolved disks mean that we cannot make a reliable assump- 
tion as to the composition of the grains in an unresolved 
disk. In addition, the poor fit to e Eri's SED shows that the 
model may not in general contain all of the physics necessary 
to account for the observations. Our estimates of the disk 
size are based on assumption that the dust grains are solid 
(as for Fomalhaut), but if the grain properties axe different, 
then the size estimates will be incorrect. 

Future progress in this field depends on spatially resolv- 
ing disks, as this allows a direct measure of disk size and thus 
measurements of the dust composition. Our results indicate 
that there are probably a few more targets that can be re- 



solved using SCUBA, given sufficient integration time. Also, 
ground based mid-IR observations also have the ability to 
resolve some disks. However, large numbers of resolved disks 
may have to wait until SIRTF and SOFIA become available. 



8 CONCLUSION 

We have presented new observations and modelling of Vega 
excess disks with sub-mm data. Our observations expand the 
sample of Vega excess stars detected in the sub-mm from 5 
to 10, with a further 4 objects which were not detected but 
have useful upper limits. We have fitted the observed SEDs 
with models based on realistic dust grain composition and 
size distribution. We find that dust grain composition varies 
significantly between different objects, with younger disks 
having more porous grains. For e Eri, our model fails to fit 
the data unless there are less small grains in the system than 
expected from a coUisional cascade. This discrepancy may be 
due to an inadequacy in our model of the size distribution, or 
the result of an external influence such as a planet. For the 
unresolved targets, disk size is estimated Eissuming solid dust 
grains, but given the observed diversity in dust composition 
and the problems fitting e Eri, these size estimates may not 
be reliable. 

The James Clerk Maxwell Telescope is operated by 
the Joint Astronomy Centre on behalf of the UK Particle 
Physics and Astronomy Research Council. We wish to thank 
the staff at the JCMT for enabling these observations, and 
an anonymous referee for their helpful and positive com- 
ments. 

This paper has been typeset from a T^jX/ file prepared 

by the author. 
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